Optimizing number squeezing when splitting a mesoscopic condensate 
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We optimize number squeezing when splitting a mesoscopic Bose Einstein condensate. Applying 
optimal control theory to a realistic description of the condensate allowed us to identify a form of the 
splitting ramp which drastically outperforms the adiabatic splitting. The results can be interpreted 
in terms of a generic two-mode model mapped onto a parametric harmonic oscillator. This optimal 
route to squeezing paves the way to a much longer phase coherence and atom interferometry close 
to the Heisenberg limit. 
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Confined atom interferometers using Bose-Einstein 
condensates (BECs) offer new prospects for matter wave 
interferometry [1 and precision measurements. Optical 
dipole traps [2], atom chips [3], and radio-frequency (rf) 
potentials [4] [5] provide powerful tools which enable co- 
herent manipulation and interference as demonstrated in 
a series of recent experiments [6j [7] [8] [10] [11] . 

Atom interferometers based on BECs usually suffer 
from the nonlinearity originating from atom-atom inter- 
actions, which leads to phase diffusion [12]. A possi- 
ble way out is to seek for narrow number distributions 
of the split condensate- i.e., squeezed states, which are 
very powerful in precision measurements [13] [HI EI HE] • 
This can be achieved by adiabatic splitting, where the 
nonlinear interaction favors narrow number distributions 
[8] [10] [17] . The disadvantage of this scheme is the long 
time needed for the splitting process, within which tech- 
nical noise and additional phase diffusion might threaten 
the interferometer performance [18] [19] [20] . 

In this paper, we show that splitting protocols, based 
on optimal control theory (OCT), allow efficient num- 
ber squeezing on a much shorter time scale and dras- 
tically outperform adiabatic splitting. We first inves- 
tigate the OCT problem in the framework of a simple 
two-mode model, leading us to an intuitive interpreta- 
tion of the control strategy. The predictions of the sim- 
ple model are verified for a realistic experimental setting 
by applying OCT to the many-body problem, within 
the framework of multiconfigurational time-dependent 
Hartree equations for bosons [MCTDHB(2)] [21 . The 
fringe visibility is significantly enhanced in case of opti- 
mized splitting, which renders this scheme ideal for atom 
interferometry. 

Splitting a Bose-Einstein condensate is achieved by 
changing the confinement potential smoothly from a sin- 
gle well to a double well, as schematically shown in the 
inset of Fig. 1. We will assume that the condensate wave 
function is modified only along a single spatial direction 
x. To describe properly the fragmentation of the BEC 
into two spatially separated condensates, we need at least 
two wave functions </)l,r(x), which we will refer to as 
left and right orbitals, together with an additional part 
that describes how the atoms are distributed among these 



two orbit als. Close to the splitting point, where the two 
orbitals become spatially separated, the system can be 
approximately described by a generic two-mode model, 
characterized by the Hamiltonian [TTJ [22] 

H=-^E j J x + 1 -E c J 2 z . (1) 

Here, the Josephson energy Ej, which is proportional to 
the energy overlap of the orbitals, accounts for tunneling 
whereas the charging energy E c accounts for the non- 
linear coupling of the atoms. J x and J z are pseudospin 
operators associated with these couplings [TT] [22s J x 
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FIG. 1: (Color online) Time evolution of atom number fluc- 
tuations for Ej(t) — Ej(0) exp[— ^] and for different decay 
constants r, as computed within the generic two-mode model 
[T71 [22] . With increasing r, the splitting process is more 
adiabatic and the fluctuations in the final state are lowered. 
The symbols report results of our OCT calculations for dif- 
ferent numbers of atoms and splitting times T, and show that 
OCT can significantly outperform the more intuitive quasi- 
adiabatic scheme for the exponential turning off. The square 
(diamond) symbols correspond to solutions obtained within 
MCTDHB for N = 100 (N = 500) (see Fig. 3). In case of 
MCTDHB we take the E c value after splitting for compar- 
ison. Inset: Schematics of the BEC splitting process. By 
ramping up a double-well potential, the initial ground-state 
wave function of a single well becomes nonadiabatically split 
into two parts, denoted as left and right orbitals. 
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promotes an atom from the left to the right well, or vice 
versa, and J z measures the atom number difference be- 
tween the two wells. 

In experiments, Ej is controlled indirectly by variation 
of the confinement potential. Here we use Ej itself as a 
control parameter in order to grasp the essential features 
of the control strategy, and will lift this assumption later. 

When tunneling dominates over the nonlinear interac- 
tion, Ej ^> E C N, all atoms reside in the bonding orbital 
(j) g = ^(</>l + resulting in a binomial atom number 
distribution with fluctuation Ano- When the tunnel cou- 
pling is reduced, the nonlinear coupling favours localiza- 
tion of the atoms in one of the wells. The state of lowest 
energy is a superposition of different atom number states 
with smaller than binomial number fluctuations. In the 
limit of very small tunnel coupling, Ej — > 0, the ground 
state is such that half of the atoms reside in the left well 
and the other half in the right well, and there are no atom 
number fluctuations, An = 0. 

To split the condensate, one starts from a state with 
Ej ^> E C N and then turns off the tunnel coupling. A 
quasiadiabatic splitting corresponds to an exponential 
decrease of Ej, as shown in Fig. 1 for different decay 
constants r and 100-1000 atoms, corresponding to real- 
istic experimental conditions. One observes that for a 
slowly varying Ej the system evolves almost adiabati- 
cally and finally ends up in a state with small atom num- 
ber fluctuations. For faster splitting, the system can no 
longer follow adiabatically and becomes frozen in a state 
with substantially larger number fluctuations (less num- 
ber squeezing). By increasing r by a factor of 10, the 
number fluctuations in the final state drop by a factor of 
approximately 2. Thus, efficient atom number squeezing 
comes at the price of very slow splitting. 

To improve this we are seeking for an optimal time 
variation of Ej(t) that brings the system to a number 
squeezed state with reduced atom number fluctuations 
in much shorter time. To this end, we employ OCT 
[23j [24] [25] with the goal to minimize the atom number 
fluctuations (An) 2 
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FIG. 2: (Color online) Results of OCT calculations for the 
generic two-mode model and for an atom number N — 100. In 
(upper panel), we plot the tunneling control Ej(t), in (middle 
panel) the time evolution of the atom number fluctuations, 
and in (lower panel) a density plot of the absolute square 
modulus of the atom number wave function. The bright lines 
show estimates based on the parametric oscillator of Eq. Q. 
The time interval has been chosen somewhat larger as in Fig. 3 
to make the oscillating control mechanism more visible. 



of our squeezing optimization for N=100. OCT comes up 
with an oscillating tunnel control, which leads to a drastic 
reduction of the number fluctuations in comparison with 
the exponential turning-off. 

We next investigate why the oscillating OCT tunnel 
coupling drastically outperforms the more intuitive expo- 
nential decay. We show how to qualitatively understand 
the mechanism of the control. 

For large atom numbers iV, the time evolution of the 
generic two-mode model can be approximately described 
by a harmonic oscillator [T7] 
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in the state at the final time T. Within the framework of 
OCT, J is called the cost function, which is minimized 
under the constraint that the system's time evolution is 
governed by the Schrodinger equation. This is done by 
using Lagrange multipliers to turn the constrained mini- 
mization problem into an unconstrained one, as discussed 
in some length in Refs. [24] [25]. 

Our OCT calculations within the two-mode model are 
summarized by the circular (for N=100) and triangular 
(for N=1000) symbols in Fig. 1. One clearly sees that 
the optimization of number squeezing works within a 
wide range of splitting times, and Ej(t) sequences found 
in OCT perform approximately one order of magnitude 
faster in comparison to the standard exponential one. 

Fig. 2 shows details to one of the Ej (t) time sequences 



where C(k) is the atom- number wave function and k is 
the number difference between the left and right well, 
which is treated as a continuous variable. Introducing 
the annihilation and creation operators a and a) for the 
harmonic oscillator [26] . we can cast the Hamiltonian of 
Eq. (|3| in the form 



H = 2Ej/N 



a) a + 



1 



E C N 2 
16Ei 



(a 2 +a^ 2 ) 



• (4) 



Here Ej/N = Ej/N + E C N/S is the renormalized oscil- 
lator frequency. Equation Q is the Hamiltonian for a 
parametric harmonic oscillator [27] with resonance fre- 
quency UJ res « 2Ej /N + E c N/4. 

When the oscillator is driven with approximately twice 
the resonance frequency, which is approximately the pe- 
riod of the OCT control in Fig. 2, the width of the initial 
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ground state wavepacket starts to oscillate and becomes 
strongly squeezed. The predictions of the parametric os- 
cillator model are plotted in Fig. 2 for the oscillating 
control field indicated in the upper panel: indeed, as can 
be observed from the number fluctuations in the middle 
panel, the envelope of An(t) decreases in a fashion similar 
to the results of the OCT calculations. Thus, if we turn 
off Ej at the lower turning point of An(T), we freeze the 
system in a squeezed number state [although the detailed 
freezing sequence of Ej(t) happens to be more complex 
in case of OCT]. 

We now address the question whether our findings 
would prevail in case of a more realistic modelling of the 
many-body splitting process. Contrary to the two-mode 
model, where all details of the condensate orbitals are 
embodied in Ej and E c , in general the control of Ej and 
E c is indirect through the condensate orbitals, which, 
in turn, can be manipulated by means of the confine- 
ment potential V\(x). Here, X(t) is a control parameter 
that describes the variation of the confinement potential 
when changing the external parameters, such as currents 
through the microtrap wires or frequency and strength 
of additional rf fields [4 . 

For a realistic description of the splitting process we 
employ MCTDHB(2) [21 , where the orbitals are deter- 
mined self-consistently from a variational principle. This 
approach accounts in a natural manner for both the atom 
number fluctuations and the orbital dynamics. We raise 
the issue, if the nonadiabatic dynamics of the orbitals al- 
lows us to exploit the same control mechanism as before. 

In our calculations we use parameters typical for cigar- 
shaped potentials with a few 100-500 atoms, split along 
the transverse direction, similar to recent atom chip [6j[7] 
or squeezing experiments [10 (see caption of Fig. 3 for 
details). We first choose X(t) such that the tunneling de- 
cays approximately monoexponentially. Figure 3 shows 
that for such control the number fluctuations decay in a 
fashion similar to the generic two-mode model and the 
number squeezing in the final state is rather low. 

Squeezing optimization is again performed within the 
framework of OCT [28]. Fig. 3 shows the details of our 
OCT calculations. The density plot of the orbitals clearly 
shows that, as in the two- mode model calculation, the 
condensate is first brought to oscillations within the two 
wells, resulting in an oscillating tunnel coupling. In this 
regime the atom number fluctuations first oscillate wildly 
and then significantly drop. These results are in good 
agreement with the generic two-mode model; see Fig. 1. 

To turn off the condensate oscillations after the squeez- 
ing optimization, we introduced an additional optimiza- 
tion step for the trapping of the orbitals, similarily to our 
previous work on the optimization of the Gross-Pitaevskii 
equation [25 . With this the orbitals are brought to an 
almost complete halt, as evidenced by the stationary evo- 
lution at later times [29] . 

We obtain similar results for larger atom numbers (500 
atoms) and different splitting times. This makes us be- 
lieve that the oscillating control for achieving a high de- 
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FIG. 3: (Color online) Optimized versus unoptimized split- 
ting, as computed within MCTDHB(2). For the trap we use 
splitting by rf potentials 4 , with cut = (2n)2 kHz transversal 
frequency, and the control field A relates to the amplitude of 
the rf field B Y { = (0.5 + 0.3 A) G. For an interaction strength 
[21] (chemical potential) of g ~ 1.16 Hz (/x ^ 2 kHz for the 
unsplit trap), OCT improves squeezing considerably as com- 
pared to the exponential case. In the calculations we use 
gerade and ungerade orbitals (/> g , e , rather than as they 

allow to fully exploit the symmetry of the confinement po- 
tential. The scheme is quite robust to typical experimental 
noise of 30 /xG (black lines, indistinguishable from OCT re- 
sults). Panel (a) reports the control fields and atom number 
fluctuations for N = 100 for the exponential (dashed lines) 
and optimal control (solid lines), and panels (b) and (c) the 
time evolution of the gerade and ungerade orbitals. 
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FIG. 4: (Color online) Useful squeezing £r for N = 100 (dark 
lines) and N = 500 (bright lines) for the exponential (dashed 
lines) and optimal control (solid lines). The optimized solu- 
tions stay well beyond the standard quantum limit (£r = 1), 
in sharp contrast to the unoptimized cases. The lines in the 
lower part indicate the Heisenberg limits = y/2/N for the 
investigated TV's. 



gree of squeezing is of general nature and that the simple 
two-mode model provides a proper description for the 
underlying physics. 

In the free-time evolution after splitting, an atom num- 
ber superposition state with finite An undergoes a spread 
of evolution rates due to the nonlinear atom-atom inter- 
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actions [T2] . This "diffusion" of the relative phase with 
time degrades the coherence a := 2Re{a^ R aL) /N of the 
condensates, directly observable as the fringe contrast 
[30] . Here, cll (clr) represents the mode operator for the 
left (right) condensate. To quantitatively analyze the 
improvement of OCT, we calculated the coherence 30 ms 
after splitting for N — 100 (500) atoms, and found values 
of above 80% (90%) for OCT, to be contrasted with the 
values of close to zero (50%) for exponential splitting. 
Thus our squeezing protocol strongly improves the phase 
coherence for a long time after splitting. 

Our OCT protocols will also be useful in experi- 
ments based upon the measurement of an atom num- 
ber difference between the wells. The phase sensitiv- 
ity is then quantified by the factor of useful squeezing 
£r = 2An/\^Na [14]. It displays that squeezing en- 
hances the sensitivity below the shot-noise limit £r = 1 
and is (Heisenberg) limited from below by y/2/N [T3] . 



Furthermore, £r < 1 is a sufficient criterion for the pres- 
ence of entanglement between the N atoms, signifying 
a type of entanglement which is useful as a resource 
[TO] I3T] . It is apparent from Fig. 4 that OCT squeez- 
ing achieves a phase sensitivity close to this fundamental 
limit of quantum measurements as well as the creation of 
N-body nonseparable states. 

In conclusion, we have demonstrated that non- 
adiabatic condensate splitting following OCT allows for 
a very efficient number squeezing on short time scales 
leading to a strongly enhanced phase coherence, thus ren- 
dering the technique powerful for interferometry applica- 
tions. 
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